Sharp bends in photonic crystal waveguides as nonlinear Fano resonators 



in 
o 
o 

(N 



Andrey E. Miroshnichenko and Yuri S. Kivshar 
Nonlinear Physics Centre and Centre for Ultra-high bandwidth Devices for Optical Systems ( CUDOS), 
Research School of Physical Sciences and Engineering, 
Australian National University, Canberra ACT 0200, Australia 

We demonstrate that high transmission through sharp bends in photonic crystal waveguides can 
be described by a simple model of the Fano resonance where the waveguide bend plays a role of 
a specific localized defect. We derive effective discrete equations for two types of the waveguide 
bends in two-dimensional photonic crystals and obtain exact analytical solutions for the resonant 
transmission and reflection. This approach allows us to get a deeper insight into the physics of 
resonant transmission, and it is also useful for the study and design of power-dependent transmission 
through the waveguide bends with embedded nonlinear defects. 

I. INTRODUCTION 

o . 

^ Photonic crystals (PCs) are artificial dielectric structures with a periodic modulation in the refractive index that 

create a range of forbidden frequencies called a photonic band gapi. The existence of the photonic band gap can 

change dramatically the properties of the light allowing the realization of ultra-compact and multi-functional optical 

devices. One of the most fascinating properties of photonic crystals is their ability to guide electromagnetic waves 

in the narrow waveguides created by a sequence of line defects, including light propagation through extremely sharp 

waveguide bends with nearly perfect power transmission^. It is believed that the low-power reflection of sharp 

waveguide bends in photonic crystals is one of the most promising approaches to combine several devices inside a 

I , . compact nanoscale optical chip. 

The main advantage in achieving low radiation losses for the light transmission through sharp waveguide bends is 

based on the existence of the photonic band gap allowing to confine light inside a narrow defect waveguide due to the 

effect of the resonant Bragg scattering, reducing only the reflection losses of the bend design. Recent studies addressed 

the issue of an improved design of the sharp waveguide bends in two-dimensional photonic crystals and suggested 

that the transmission losses can be less then 5%2ii£i|jiiL£ii!ii2iiiii2iIii. However, most of those studies are based on 
t- H ... . 

direct numerical simulations and neither provide a deep understanding of the physics of the enhanced transmission of 

the sharp waveguide bends nor suggest simple approaches to design the waveguide bends with required transmission 
' properties. 

t-H , The main purpose of this paper is twofold. First, we demonstrate that transmission of electromagnetic waves 
C*~) • through sharp bends of photonic crystal waveguides can be described by a simple discrete model where the waveguide 
bend plays a role of a specific localized defect. By employing the semi- analytical approach based on the Green's 
function formalismi2ii£ii& ; we derive effective discrete equations which allow us to describe the bend transmission 
as a special case of the so-called Fano resonance^, recently discussed for the wave propagation in discrete chains^. 
We demonstrate how to introduce the effective discrete model for the Fano resonance by selecting two common 
designs of the waveguide bends in two-dimensional photonic crystals created by a lattice of rods, and obtain exact 
analytical solutions for the resonant transmission and reflection. We show that asymmetric shapes of the transmission 
curves observed for the waveguide bends can be understood in terms of the Fano resonance which originates from 
the interaction between continuum waves and an effective localized state associated with the waveguide bend that 
provides an additional propagation channel for the wave transmission and, therefore, leads to the constructive or 
destructive interference. Second, we show that this approach allows us not only get a deeper physical insight but it 
is also useful to study the nonlinear transmission through the waveguide bends with nonlinear defects. 

The paper is organized as follow. In Section 2 we derive the effective discrete equations for the case of two- 
dimensional photonic crystal waveguides based on the Green's function technique. We employ the approach earlier 
suggested by Mingaleev et but take into account the effect of the long-range interaction only in the vicinity 

of the waveguide bend. This allows us to derive an effective discrete Fano model where the waveguide bend plays 
a role of a special localized defect, and also suggest an effective way for predicting and controlling the properties of 
different types of the waveguide bends. Using this model, in Section 3 and Section 4 we study the wave transmission 
in the linear and nonlinear regimes, respectively. Section 5 concludes the paper. 
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FIG. 1: (a) Spatial structure of the coupling coefficients J nm (uj) of the effective discrete model 10 at u) — 0.4 x 2irc/a. (b) 
Dependence of the specific coupling coefficients (marked) on the frequency u. 



II. EFFECTIVE DISCRETE MODEL 



First, we derive an effective discrete model for the wave transmission through a sharp waveguide bend by employing 
and modifying the conceptual approach suggested earlier by Mingaleev et ah 14 ^ . We consider a two-dimensional 
photonic crystal created by a square lattice (with the period a) of dielectric rods in air (eb g = 1). The rods have the 
radius r = 0.18a and the dielectric constant e ro d = 11.56. We study in-plane light propagation in this photonic lattice 
described by the electric field E(x, t) = exp(— iu)t)E(x\uS) polarized parallel to the rods, and reduce the Maxwell's 
equations to the scalar eigenvalue problem 

V 2 + n 2 e(x) 



E{x\uj) = . (1) 

For given parameters, this square lattice of rods is known to possess a large TM band gap (38%) between the 
frequencies cu = 0.303 x 2itc/a and u> — 0.444 x 2irc/a. 

We create a waveguide by replacing some of the lattice rods by the defect rods with the radius r<2, or simply by 
removing some rods of the lattice. To describe the structure with defects, we decompose the permittivity function 
e(x) into a sum of the periodic part and the defect-induced contribution, e(x) = e p (x) + Se(x), and rewrite Eq. Q in 
the integral formi^, 

E(x\u>) = {^) 2 J d 2 yG(x,y\Lo)5e(y)E(y\Lo) , (2) 

where G(x, y|w) is the standard Green's function. If the radius of the defect rod r& is sufficiently small, the electric 
field E(x\u>) inside the rod is almost constant, and the integral can be easily evaluated. This allows us to derive 
a set of discrete equations for the electric field 

E n ,m = 2, Jn-k,m-l(u)Sek,lEk,l (3) 



k,l 



where 



Jn,m(w) = J d 2 yG(x„,x m +y\w) (4) 

rd 

are the frequency-dependent effective coupling coefficients, 

&€n,7n — £n,m ^rod ) (5) 
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FIG. 2: Schematic view of two designs of the waveguide bends studied in the paper. Empty circles correspond to removed 
rods, dashed lines denote the effective coupling. Yellow circles mark the defects with different dielectric constant ed(|i?|), which 
can be nonlinear. 



is the defect-induced change of the lattice dielectric function, where e n m is the dielectric constant of the defect rod 
located at the site (n,m). 

In general, the effective coupling coefficients |J n ,m(w)| decay slow in space, as shown in Fig.^a). This slow decay 
introduces effective long-range interaction (LRI) from site to site of the waveguide, which becomes crucially important 
for the waveguide bends. In reality, we define a finite distance L of this (formally infinite) interaction by assuming 
that all coupling coefficients with the numbers |n — k\ > L and \m — 1\ > L vanish. As demonstrated in Ref>i£, the case 
L = 6 gives already an excellent agreement with the results of the finite-difference time-domain numerical simulations. 

Unlike the previous studies of the effective discrete model, here we are interested in the light propagation through 
two types of the 90° waveguide bends shown in Fig. [3 created by removed rods and additional defect rods with 
dielectric constant td placed at the corner; the defect rod can also be nonlinear. Due to the effective long-range 
interaction, the corner in the bend waveguide can be considered as a special type of a defect; this resemblance is 
indeed correct and it explains the Fano-type resonances in the transmission as discussed below. 



III. LINEAR TRANSMISSION 



In this section, we consider the linear transmission through the waveguide bends when all corner defects are assumed 
to be linear. In order to show analytically that the waveguide bend generates an effective defect state and it can be 
described by a discrete model of the Fano resonance, we proceed in two steps making simplifications in the general 
model (01. For simplicity, we consider the waveguide bend of the type A, as shown in Fig. [2 

As the first step, we reduce the effective length of the long-range interaction in the straight parts of the waveguide 
bend taking it as L = 1. This approximation corresponds to the familiar tight-binding approximation that takes 
into account only the coupling between the nearest neighbors in a discrete chain. However, we take into account the 
nonlocal coupling near the sharp bend, and this requires to introduce an additional coupling of two sites around the 
corner as shown in Fig. [21 for the bend of the type A. It turns out that this simple approximation allows us to describe 
qualitatively all important properties and the key physics of the waveguide bend. 

To derive the effective discrete model corresponding to the bend design shown in Fig. El (left), we remove two 
semi-infinite rows of rods and write the dielectric function in the form, 

£n,m = £bg ! for n < — 1 and m = or n = and m > 1 . (6) 

At the corner of the waveguide bend, we place a defect rod with the dielectric constant eo,o = Equations |J5J for 
the electric field inside the defect rods can be written explicitly as 

(1 — Jo,o8co)E n fi — 5eoJi,o(E n +i,o + E n _ifi), n < — 1, 

(1 — Jofifi e o)Ea,m — <5eo^O,l(£'0,m+l + -E-m,n-l), m > 1, 
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(1 - Jo,o5e )-B_x,o = Ji,o(SeiE ,o + (5e ^- 2 ,o) + Ji, iSe E 0tl , (7) 
(1 — Jo,oSei)E ,o = 5e (Jo.iEo,i + Ji^-B-i.oji 
(1 — Jo,oSeo)Eo,i = Jo,i(SeiE 0i o + 5eoE ^) + Ji : iSe E^ifi, 

where Seg = e^ g — e ro( j and 5e\ = e d — e ro d- Importantly, the resulting set of coupled equations present a discrete 
model that can be compared with the discrete models earlier studied in Refii&, for which the existence of the Fano 
resonance has been demonstrated analytically. 

First two equations in the system (JZJ) allow us to obtain the dispersion relation for the bend waveguide far away from 
the corner. Due to the symmetry of the photonic crystal, the coupling terms Jo.i and Jifl coincide and, therefore, we 
obtain the waveguide dispersion in the form 

1 — 5e Jo,o fo s 

cosfc = — — , (8) 

oe J ,i 

where k is the wavenumber for the waves propagating along the waveguide. Other equations in the system J7J allow 
us to calculate the transmission coefficient of the waveguide bend, 

_ 4a 2 sin 2 fc 

\b(c-b)\*' W 

where we use the notations 

a = (Ju + Jo,i5ei - JoMJi,iSei)J ,i5el , b = ( J ,o + exp(ifc) J ,i - Ji,i)Seo - 1 , 

c = ( J o,o - 2J o,i + Jo.oJi.i + exp(ifc) J .q J .i)5e Q 6ei - Jo.qSci . (10) 

Important information about the resonant transmission follows from the study of zeros of the transmission coefficient 
@. Except for the band edges (k = 0,7r), the transmission coefficient vanishes when a = 0, i.e. when 

Jo,oJi,iSei = Ju + Jo^Sei. (11) 

As the second step, we approximate the dependence of the coupling coefficients J„.m in the effective discrete model 
on the frequency u>, considering the specific results presented in Fig. db). In the frequency interval [0.32,0.42], we 
find the following approximation 

^o,o = ijJ - ud, Jo,i = C, Ji,i = V Jo,o = V{uj — ujd), (12) 

where ujd is the frequency of a single embedded defect (in our case, u>d ~ 0.4), C and V are constants (C w 0.035 and 
V f» 0.07). By substituting Eq. 112[l into Eq. (| 1 1 p . we obtain the quadratic equation 

8eiVoj 2 - {25tiu d + 1)Vlj + (VSenvj + Vto d - C 2 8ei) = 0, (13) 

which has two solutions 



UJ F =UJ d + ± 
L0t\ 



AC 2 



V 
464 



1/2 

(14) 



According to these results, there exists a possibility for two perfect reflections and, therefore, two Fano resonances. 
This is due to the linear frequency dependence of two coupling terms <|12[) . the on-site term Jo,o and the cross-coupling 
term J\^\. In our model, we have only one free parameter — the dielectric constant e d of the defect rod at the corner. 
Then, when 8e\ ~ at list one zero of the transmission lies inside the propagation spectrum JSJ, loq < up < uj v . 
As a result, we predict that when e d ~ e ro d the perfect reflection through the waveguide bend of the type A should 
be observed, and numerical results confirm this prediction (see Fig. |3J). We plot the transmission coefficient for two 
different values of the interaction, L = 1 and L = 6. This shows that the tight-binding approximation (L = 1) still 
gives reasonable results in average, and it works especially well near the Fano resonance where our model works nicely. 

Both the effective model and calculation of the transmission coefficient for the waveguide bend of the type B (see 
Fig|2J are similar. The only difference between these two types of the bend design is the existence of additional 
coupling terms to the the defect state. As was shown in Ref^, by increasing the number of the coupling terms we 
only shift the position of the Fano resonance, either transmission or reflection. In our case, the coupling is small 
(V <C 1), so that the renormalization is negligible. For the waveguide bend of the type A, the Fano resonance 
manifests itself as the perfect reflection only, whereas the perfect transmission lies outside the waveguide spectrum. 




FIG. 3: Transmission coefficient of the waveguide bend for different vaiues of the dielectric constant ed- The Fano resonance 
is observed when the value of the dielectric constant of the defect rod td approaches the value of the dielectric constant of the 
lattice rod e ro( j. The plot = tbg corresponds to the case when a rod is removed from the bend corner, whereas the plot 
£d = e ro d corresponds to the case when the lattice rod remains at the corner. 



By analyzing the transmission coefficient for the type B bend, we can show that in this case the perfect transmission 
may move to the waveguide spectrum when ed ~ £bg- This result coincides with the well-known result of the perfect 
transmission through the waveguide bend of the type B when the rod is removed at the corner, ed = eb^^- 

According to Fig. [31 by increasing the dielectric constant of the defect rod at the corner, we can achieve the 
perfect resonant reflection. However, the transmission itself becomes very low. However, it was shown in Ref^& that 
when the coupling to the defect is very small the frequency of the perfect transmission is located very close to that 
of the perfect reflection; this results in a narrow and sharp Fano resonance. We employ this idea for the waveguide 
bend of the type A and replace one lattice rod by a defect rod outside the corner, as shown in the insert of Fig. 0] 
In this case, we obtain the perfect transmission for a particular frequency by varying the dielectric constant of the 
defect rod, see Fig. 0] Here, there exist two Fano resonances. One of them is broad, and it manifests itself as the 
perfect reflection only, being similar to the case of Fig. Etc) ; it can be treated as the background transmission. The 
second resonance possesses a sharp asymmetric profile with both resonant transmission and reflection. Our analysis 
shows that for this design the so-called double Fano resonance can exists due to the specific frequency dependence 
of the coupling coefficients. By varying the dielectric constant of the defect rod placed outside the corner, we shift 
the asymmetric resonance whereas preserving the other one. When we simply remove the rod outside the corner 
(i.e. at ed = fibg), these two resonances are located very close and interact with each other. As a result, the perfect 
transmission is accompanied by two perfect reflections, and this effect can be use to design a very efficient filter based 
on the bend waveguide transmission. 
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FIG. 4: Transmission coefficient through the waveguide bend with a (green) defect rod placed outside the corner. In this 
case, there exist two Fano resonances, one of them is characterized by an asymmetric profile and corresponds to the perfect 
transmission. 
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FIG. 5: Nonlinear transmission calculated for two types of the waveguide bends shown in Fig.|3 In both the cases, the Fano 
resonance is observed as the perfect reflection. The waveguide bend of the type B allows the perfect transmission that can be 
also tuned. 



IV. NONLINEAR TRANSMISSION 



Finally, we apply our effective discrete model to the case of the nonlinear transmission through the waveguide bend 
with embedded nonlinear defect rods. In this case, we assume that the (green) defect rods placed at the bend corner 
(see Fig. EJ) possesses a Kerr-type nonlinearity 

e d (\E\) =e d + \\E\ 2 . (15) 

For definiteness, in our numerical simulations we take the value td = 1-96 and use the rescaled coefficient A — 1, these 
data should correspond to some polymer materials. In the nonlinear regime, the transmission of the waveguide bend 
depends on the intensity of the incoming light. This gives us an additional possibility to control the transmission 
properties of the waveguide bend by changing the properties of the nonlinear defect. 
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Our analysis shows that the presence of nonlinear defect does not remove the Fano resonance itself, but instead 
it shifts the position of the resonant frequency. As a result, we can tune the value of the resonant frequency by the 
input light intensity, the resonant scattering can be observed for almost all frequencies of the waveguide spectrum^. 
Figure [S] shows the power-dependent transmission for two types of the waveguide bends with embedded nonlinear 
defects at the corner. For the type B bend, the perfect transmission is observed as well. The intensity-dependent 
transmission allows us to control the light propagation through the waveguide bend from 0% up to 100% by simply 
tuning the light intensity. From another hand, we can achieve 100% transmission thought the type B bend for almost 
all frequencies from the spectral range by choosing a proper intensity of the incoming light. 



V. CONCLUSIONS 



We have analyzed the conditions for high transmission through sharp waveguide bends using the effective discrete 
equations derived for two-dimensional photonic crystal waveguides. We have demonstrated that the physics of this 
effect can be understood with the help of a simple discrete model of the Fano resonance where the waveguide bend 
plays a role of a specific localized defect. Using this model, we have obtained exact analytical solutions for the resonant 
transmission and reflection of two types of the waveguide bends, in both linear and nonlinear regimes. We believe 
our approach would allow to get a deeper insight into the physics of resonant transmission through waveguide bends, 
and it can be useful for understanding other types of resonant effects in two- and three-dimensional photonic crystal 
waveguides and circuits. 

The authors acknowledge useful discussions and collaboration with Sergei Mingaleev. 
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